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Abstract. Using generating functional and replica techniques, respectively, we study 
the dynamics and statics of a spherical Minority Game (MG), which in contrast with 
a spherical MG previously presented in yQ displays a phase with broken ergodicity and 
dependence of the macroscopic stationary state on initial conditions. The model thus 
bears more similarity with the original MG. Still, all order parameters including the 
volatility can computed in the ergodic phases without making any approximations. 
We also study the effects of market impact correction on the phase diagram. Finally 
we discuss a continuous-time version of the model as well as the differences between 
on-line and batch update rules. Our analytical results are confirmed convincingly by 
comparison with numerical simulations. In an appendix we extend the analysis of pQ 
to a model with general time-step, and compare the dynamics and statics of the two 
spherical models. 



PACS numbers: 02.50.Le, 87.23.Ge, 05.70.Ln, 64.60.Ht 
E-mail: galla@ictp.trieste.it, sherr@thphys.ox.ac.uk 

1. Introduction 

The introduction of continuous versions of spin systems with discrete degrees of freedom 
is common in statistical mechanics and has led to a better analytical understanding 
of many basic features displayed by various spin models. The appeal of continuous 
models lies in the fact that calculations are simplified considerably compared to their 
counterparts with discrete Ising spins, so that in many cases a full analytical solution 
of the models with continuous spins becomes possible. The general procedure is 
to replace the discrete Ising-like variables G { — 1,1} of a given spin system by 
continuous variables (pi G 1R. If we assume that the system is composed of N such spins 
(i = 1,...,N), then often a so-called spherical constraint is imposed on all possible 
configurations contributing to the partition function: <fif = N. While the partition 
function of a system with Ising spins is a sum over configurations corresponding to 
the corners of an iV-dimensional hypercube, this sum is effectively replaced by an 
integral over the surface of the sphere running through these corners. The hope is 
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that the simplified, spherical model still captures the relevant physical features of the 
original model and allows to shed some light on its basic phenomenology despite the 
simplification. 

The first example of a spherical model is the one of a ferromagnet introduced more 
than fifty years ago by Berlin and Kac It is the continuous analogue of the Ising 
model, and can be solved exactly in all dimensions, and exhibits a phase transition 
corresponding to the onset of spontaneous magnetization. Spherical counterparts of the 
Sherrington-Kirkpatrick model have played and are playing an important role in the 
study of spin-glasses, and similar mean-field models with higher-order couplings between 
the continuous spins are able to capture some of the essential features of structural 
glasses. Recently a spherical version of the Hopfield model of neural networks has been 
considered in jHj, and a spherical Minority Game (MG) was introduced in pQ. The 
present work is a continuation of pQ, and we will here devise an alternative way of 
formulating a spherical limit of the MG. 

The Minority Game [I] is a simplified model of interacting traders in a financial 
market and has become one of the most studied models in the econophysics community 
|3 El Ej • From the point of view of statistical mechanics the MG is a fully connected 
mean-field disordered system. Analytical progress can therefore be made using the 
standard tools of equilibrium and non-equilibrium statistical mechanics, such as replica 
theory and generating functionals. With these methods an ergodic/non-ergodic phase 
transition has been identified in the MG, and static order parameters in the ergodic 
regime computed within a replica-symmetric ansatz agree with numerical simulations. 
The dynamical approach furthermore allows one to formulate a self-consistent theory 
for all regimes of the game in terms of a single effective-agent process jH] . This effective 
problem is fully equivalent to the coupled multi-agent dynamics of the MG in the 
thermodynamic limit and analytical results for persistent order parameters in the ergodic 
phase can be obtained from this effective process in agreement with the results derived 
from the statics. These solutions are now regarded as exact. A proper understanding of 
the non-ergodic phase is however still awaited. Moreover, for the conventional MG the 
so-called 'market volatility' (one of the MG's main observables) cannot be computed 
exactly and results are here restricted to approximations even in the ergodic phase. 

In [l j a spherical version of the MG was presented, offering the possibility of 
a more complete solution. The non-linear dynamics of the conventional MG was 
replaced by effectively linear equations, resulting in an exactly solvable harmonic model. 
Despite the simplicity of the update rules this model was shown to exhibit a non-trivial 
phase diagram with three distinct phases and different types of phase transitions. All 
stationary observables (including the volatility) can here be computed in all phases 
without making any approximations. The model discussed in 1] however has only 
limited similarity compared to the original MG. For example the behaviour of the 
volatility turns out to be quite different in both models, and crucially the spherical 
model presented in pQ does not exhibit any phases in which the macroscopic dynamics 
and the resulting stationary states depend on initial conditions. This is in sharp contrast 
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with the original MG, in which initial conditions have been found to be relevant in the 
non-ergodic phase and where the resulting memory effects have to a great part been 
responsible for the interest of statistical physicists in the MG. 

In this paper we will therefore consider a spherical limit of the MG, which is different 
from the route taken in pQ , leading to a model with a phenomenology closer to that of the 
original MG. In particular the model studied in this paper displays non-ergodicity and 
memory effects. Nevertheless it is exactly solvable, and we are able to derive expressions 
for the volatility in the ergodic regime without making any approximations. We here 
choose to concentrate on the stationary states of the model and do not address its 
transient dynamics. While we present some numerical results for the non-ergodic phase 
a more detailed analysis would here require a full solution for the transient dynamics in 
order to keep track of the effects of the state from which the dynamics is started. This 
is beyond the scope of the present paper but might be considered in later studies. 

The remainder of this paper is organized as follows: In section |21 we introduce the 
model, sectionElis concerned with the study of the dynamics, and in sectionEJwe discuss 
the resulting phase diagram. The statics are considered in sectional Finally we present 
a comparison of the batch model with its continuous time limit (section |HJ) and with 
its on-line counterpart (section EJ), before we draw conclusions in the final section. The 
appendix provides a discussion of some aspects of the earlier spherical MG not addressed 
in pp as well as a comparison of the two spherical models. 

2. Model Definitions 

Before turning to spherical versions of the MG, it will be useful to briefly summarise 
the definition of the MG in its conventional setup. The MG describes an ensemble of 
N agents, whom we shall label by Roman indices. At each round t of the game, each 
agent i takes a binary trading decision b^t) G {—1,1} in response to the observation 
of a publicly available piece of information I^t)- The aim of each agent is then to 
make the opposite choice to the majority of agents. In this paper we will assume 
that fi(t) is chosen randomly and independently from a flat distribution over the set 
fi(t) G {1, . . . , aN} at each time-step. This is the so-called MG with random external 
information. The key control parameter of the model is the ratio a of the number of 
possible values of the information over the number of players. We will only consider the 
case in which a does not scale with N, i.e. a = O(N ). The (re-scaled) total market 
at time-step t is then defined as A(t) = N^ 1 ^ 2 ^ bi(t). To take trading decisions each 
agent i has two strategies H ia = (R} a , ■ ■ ■ , RfJ*) G { — 1, l} aN at his or her disposal, with 
a = ±1|. These effectively act as look-up tables mapping the values \x of the external 
information onto a trading decision Rf a . If agent % decides to use strategy Si(t) G {—1,1} 
in round t his or her bid at this stage will be R^jfLy In this paper we will assume that 
all strategies are drawn randomly and independently before the start of the game and 

| The MG problem can be generalised to include more strategies per agent, but here we concentrate 
on this technically simpler but still representative case. 
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that the components {-R^} take the values ±1 with equal probability; the strategies 
then remain fixed, they represent the quenched disorder of the problem. In order to 
decide which of their two strategies to use each agent keeps track of the performance 
of both of his or her strategies by assigning virtual scores Pi a (t) based on what would 
have happened had he or she always played that particular strategy These scores are 
updated according to 



Pia(t + 1) = Pia(t) 



A (t)-^=(R^l-R^ 

Note that the minus sign in front of the square bracket ensures that strategies which 
would have produced a minority decision are rewarded. The term proportional to k 
takes into account a correction of the impact of trader i's action on the total market 
bit, and can be seen as an analogue of an Onsager reaction term; see jHJ [§] for details. 
The additional model parameter k can take values < k < 1, with k = corresponding 
to the absence of impact-correction and k = 1 to so-called 'sophisticated agents' or 'full 
impact correction'§. In the conventional MG each player then uses the strategy with 
the highest score, i.e. Sj(i) = arg ma,x a p ia (t). In order to determine which of the two 
strategies to use, it is sufficient to consider the score difference = h\pi,i(t)—Pi,-i(t)], 
then one has Si(t) = sgn and upon introduction of uji = |[Rj,i + and 
£j = IfRj^— Rj_i] the bid of player i at time t may be written as 6j(t) = CJ M ^+Sj(i)£f . 
The evolution of the {q{\ is then given by 



(1) 



»(* + 1) = »(*) - % n"« + n-^ ^e; (t) ^w - -^^sM , (2) 



where f2 = A r_1 / 2 ^ . a>j. This defines the so-called 'on-line' MG. Alternatively one may 
consider models in which the agents update the {(&(£)} only every O(N) time-steps. 
This leads to an effective average over all possible values of the external information 
and results in the so-called batch MG [Til] : 

qi(t + 1) = qi(t) - 22 JijSji 1 ) + Kasi(t). (3) 

3 

Here J tj = 2N~ 1 ^ i ■ ^ and hi = 2N~ 1 ^ i ■ Ct. 

These batch and on-line update rules provide suitable starting points for the 
introduction of spherical versions of the MG. In |JJ Eq. © was replaced by (only 
the = was considered there) 

[1 + X(t + l)] ft (t + 1) = qi (t) ~hi-J2 JiMtt) ( 4 ) 

3 

and the {A(t)} were chosen to impose a global spherical constraint ^j?i(t) 2 = Nr 2 
for all t. Although inspired by the spherical limits sometimes imposed on models of 

§ In the conventional MG the impact-correction term was found to have crucial consequences for the 
properties of the model. The game without impact correction shows an ergodic/non-ergodic phase 
transition, but no replica-symmetry breaking. In the game with impact correction (k > 0) a de 
Almeida-Thouless transition of a different type, along with the onset of long-term memory at finite 
integrated response, is found 0121 • 
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magnetic systems where Ising variables = ±1 are replaced by continuous spins fa 
subject only to the same global constraint ^ i=1 4>1 = N = J2f=i s h without the local 
constraint sf = 1, the situation in the earlier spherical MG model is at variance with 
the original MG in that the latter does not obey qf = constant. Rather it is s« = sgn[gj] 
which obeys s? = 1 for all % in the MG. In the present paper we will therefore consider the 
replacement Sj(i) — - > <f>i(t) with continuous variables {<fii(t)}, and will impose a spherical 
constraint only on the {&}, but not on the {qi}- Specifically, we will consider a batch 
model defined by 

qi(t + 1) = qi(t) ~ hi-^ Jijfyit) + «a<M*), (5) 



where 



Mt) = j@, m = [iv- 1 ^^) 2 ]"", (6) 

so that the {4>i{t)} obey the spherical constraint Yli^iW 2 = N f° r an ^ i n analogy to 
the {si(t) = ±1} in the original MG. The corresponding on-line model is defined by 

*(* + 1) = ft(o - *~ + iv-v^ef - -^ef V*w , (7) 

V L ^ v N J 

where we again impose In principle more general dependencies of the {0j} on the 
{gj} could be considered. We here chose the (pi to be proportional to the qi to guarantee 
the solvability of the resulting harmonic model. Also we assume X(t) > 0, i.e. (f>i(t) is 
taken to have the same sign as qi(t). 

An interpretation of spherical updates rules in terms of the decision making of the 
agents was given in pQ: while the agents in conventional MGs have at each time-step to 
decide which of their two strategies they wish to use, they are playing linear combinations 
of their strategies in the spherical case. The two spherical limits differ in the way this 
linear combination is determined. Moreover the models defined by Eqs. (0) and (JJJ) allow 
for runway solutions qi(t) — > ±oo (potentially resulting in an asymptotically divergent 
normalization factor A(i)), which are suppressed by the spherical constraint on the 
{qi(t)} in the earlier spherical model||. 

Let us finally, in this section, introduce one of the key observables in MGs, the 
so-called volatility a 2 . It describes the variance of the total re-scaled market bid A, and 
can be defined as the following long-time average: 

a 2 = lim rT^) 2 (8) 

T^OO ' J 

t<T 

in on-line models. In batch games an additional average over the external information 
is to be performed and one has 

a 2 = lim (aiVr)" 1 £ £ (A^t)) 2 , (9) 

T— >00 ' * ' * 

t<T [1=1 

|| Such runaway solutions, in which the |<7i(i)| increase linearly with t, are a standard feature of the 
conventional MG, corresponding to agents described as 'frozen' [HlEj- 



1/2 
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where A»(t) = W + N' 1 / 2 £\ £f<M*); for further details see also PUHUE]- Note that 
stochastic trading with the {4>i(t)} taken randomly and independently at any time step 
t from the spherical surface ^ <p 2 {t) = N would result in cr 2 = 1. This is referred to as 
the random trading limit. 

3. Generating functional analysis 

3.1. Effective single- agent process and macroscopic dynamics 

We will now proceed by studying the dynamics of the present spherical model using 
generating functionals. This approach is now standard in the context of MGs, and leads 
to a self-consistent problem for the correlation and response functions formulated in 
terms of a single effective trader equation. Due to the similarity of the present model 
with the conventional batch MG and the spherical model presented in P we can obtain 
the single effective agent equation corresponding to the process (I5I6|) by some minor 
modifications of the results of ^01 H] • One finds 

q(t + 1) = q(t) + 9{t) - a + G) a }(j)(t') + aK(j){t) + y/cen(t), (10) 



t'<t 



where <f>(t) = 

Here 6(t) is an external perturbation field introduced to generate response functions and 
f](t) is a zero- average Gaussian noise with temporal correlations characterized by the 
following covariance matrix (with D u > = 1 + Cu> for all t, t', and 1 the identity matrix): 

r* = fo(*MO>, = [(J + G)- l D(i + G T r l ] w . (11) 

The matrices C and G and the normalization factors A = {A(t)} are the dynamical 
order parameters of the problem, to be determined self-consistently upon solving 

c* = M*W)>.. G «' = ^y^))*' c * = 1 vt '*'- ( 12 ) 

The brackets (. . in (fTTlll2|) refer to averages over realizations of the process f|10)l . i.e. 
over the noise {r](t)}. The resulting self-consistent problem for {C, G, A} is equivalent 
to the original batch process in the thermodynamic limit N — ► oo. In particular, the 
physical meaning of the matrices C and G is given by 

C tv = lim N- x Y7^f)W), (13) 

i 

g -= h^i:^ (14) 

where 7TT denotes an average over the disorder, i.e. over the space of all strategy 
assignments, and (. . .) stands for an average over possibly random initial conditions 
(and/or decision noise in case stochastic trading is considered [12] )• As usual the single 
effective agent process is non-Markovian and contains coloured Gaussian noise. 

Similar to \1\ it is possible to convert the system (fTUl fTTT [TJJ) into a pair of 
explicit iterative equations for C and G and A; note that due to the non-linearity 
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in Si(t) = sgn[gj(t)] such a conversion is in general not possible for conventional MGs. 
We find 

\{t + l)C t+1 , t , - (A(t) + aK )C w = a[(l + ^-^(l + G T y l G T } tt ' 

- a[(l + G^Cjw, (15) 
X(t + l)G w - (A(t) + aK)G tt , = - a[(l + G)~ l G\ tt > + 8 W . (16) 

Eqs. (IT5| HBJ) have to be solved subject to the constraint C i4 = 1 for all t > 0. 
Furthermore one finds as in ^0] for iV — > oo that the rescaled disorder-averaged total bid 
(A(t)) is zero for all times, and that the disorder-averaged volatility is asymptotically 
given by the diagonal elements of the covariance matrix of the single-agent noise: 

a 2 = \ lim r- 1 V[(l + G)-^D{1 + G T )- 1 ] tt . (17) 

t<T 

Already at this stage one may notice that Eqs. ()15ll6j) allow in principle for non- 
1/2 



ergodic behaviour. The value A(t = 0) = liniAr^oo iV -1 ^=1 (%(0) 2 ) is n °t fixed by 
the dynamics, but is set by the initial conditions from which the dynamics is started. 
The type of stationary state reached asymptotically may hence depend on the starting 
point, similarly to what is observed in conventional MGs, and at variance with the 
behaviour of the spherical model defined by (@J), where the macroscopic dynamics at all 
times is fully fixed by the model parameter r. We will discuss the effects of different 
starts |<7i(0)| = go for all i with varying bias go below. 

3.2. Time-translation invariant ergodic states 

We will consider the system long after any initial equilibration and focus on time- 
translation invariant (TTI) solutions of the coupled equations (J151 [TH|) of the form 

lim C t+Tjt = C(t), lim G t+Tit = G(r). (18) 

t— >oo t— »oo 

In addition we make the standard assumptions of finite integrated response and weak 
long-term memory: 

X = lim Vg(t) < 00, limG M , = Vt' finite. (19) 

T<t 

The order parameter A comes out as a function of the form A(t) = Ao+Ait asymptotically 
in simulations, with Ao and Ai constants. As it turns out we will find two distinct phases 
with the above properties: 

(i) an oscillatory phase in which lim^oo \(t) remains finite, i.e. Ai = and 

(ii) a frozen phase in which X(t) grows linearly in time asymptotically, i.e. Ai > 0. 

We will refer to these phases as 'bounded' and 'unbounded' in the following, and will 
discuss them separately^. 

One may argue whether a phase with a linearly growing order parameter \(t) qualifies as time- 
translation invariant. Note however that the stationary ergodic states of the standard MG also contain 
runaway solutions of the type qi(t) — > ±00, so-called frozen agents, so that TV -1 ^ qi{t) 2 diverges in 
time. In the unbounded phase we find that \{t)/t is the stationary quantity asymptotically. 
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3.2.1. Bounded oscillatory states (Ai = 0) We here consider X(t) = Ao- Given our TTI 
ansatz (|TH|) all matrices in (fTolfTBj) become Toeplitz matrices in the stationary state, and 
hence they commute. It is then furthermore convenient to introduce Fourier transforms 
of the correlation and response functions according to 

du 



C{r) 



2vr 



-C(lu), C( W ) = ^e-T(r) 



(20) 



and similarly for G. Now upon setting A(t) = Ao in the stationary state, Eqs. (fTH| ITH jl 
translate into 

aD(uj)G(uj)* aC{uj) 



[\ (e luJ - 1) - an] C{uo) 
[\ {e iu} - 1) - <xk] G{uj) 



\l + G{u)\ 2 1 + GH 
aG{ui) 



(21) 
(22) 



1 + G(u) 

Here G{uj)* denotes the complex conjugate of G{uj). Using D(uj) = C{uj) + 2tiS(uj) as 
well as the definition of the integrated response x — J2 T ^( r ) = ^(0) we ma Y write (|2*T|) 

as 

|2 



{A (e^ - 1) - cm} |1 + G{u)\ 2 + a C{u) = 2na X S{uj). 



(23) 



From this equation it follows immediately that C(uo) can be non-zero for any w ^ 
only if e tuJ is real, i.e. for u = n. Consequently C{ui) must be of the form 
C(oj) = 2ti[cqS(uj) + ciS(u — 7r)]. Equivalently we may write 

C(r)=c + Cl (-iy. (24) 

Note that we always have cq + c\ = 1 due to the spherical constraint C(r = 0) = 1. 
Inserting this ansatz into ()21I22|) and setting u = then leads to the following relations 

l + (l-a) X 



CXKCq 



a(l + c )x . c 
+ a- 



OLKX 



[i + xY ' "i + x' i + x 

After some algebra the physically relevant solutions are identified as 

2 1 - 1 
1 



c 



1 + (1 - «)a + V4 7 ! + (1 - «0<*] 2 -4a 



2^a~ 



aV 1 + Vco - 1] 



-i 



(25) 



(26) 



(27) 



—^—r for 

a— 1 



This determines Cq and x as functions of a and k, and simplifies to Co = x 
n = 0. The expression for Co implies that a fully frozen solution (co = 1) can be realised 
only at the point given by 



a c AK) 



5 + 4k + 3 VI + 8k 
4(1 -/c) 2 ' 



(28) 



This describes a line in the (a, «)-plane, so that there can be no frozen TTI phases with 
finite integrated response and with a finite value of A = lim^oot -1 J2 T <t -M r )- We can 
therefore focus on solutions with cq < 1. These are solutions in which oscillations of 
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the correlation function persist, so that the oscillation amplitude C\ = 1 — c is strictly 
positive. Setting to = tt in (I21I22|) gives 

((_2A - aK )(l + x 'f + a )(l - c ) = 0, (29) 

(_2A - a«) X ' = 1 - -^f- r (30) 

1 + X 

Here y/ = l) r ^( r ) = measures the response of the system to persistent 

oscillatory perturbations. Looking only for oscillatory solutions (c\ > 0, Cq < 1) we 
must require (2A + «k)(1 + y/) 2 = « in this regime. Using this and (|5U|) we find the 
following solutions valid for a > a c i(/c): 

Xq = ~[a + 1 + 2y/a - an], x ' = - (31) 
2 1 + 

The persistent order parameters in the ergodic bounded oscillatory state are thus fully 
described by Eqs. (|271 l2*o1 l3*T|) . The consistency of this ansatz breaks down whenever 
X — > oo or Cq — > 1. From (J2%|) one observes that Co — > 1 at a value of a = a c i > 2, 
whereas ()27|) dictates that x remains finite for all a > 1 for any fixed value of k > 0. 
Hence we conclude that the oscillation amplitude vanishes before anomalous response 
sets in as a is lowered at a fixed value of k. We thus expect a transition to a frozen 
state at finite integrated response. 

Finally, using the fact that C(r) = cq+ (1 — co)(— 1) t the volatility in the stationary 
state can be obtained from Eq. (JT7j) as 

o 2 = i[(i + G)- 1 i>(i + G r )- 1 ] (0) = i 

Note that this is an exact result, with no approximations made at any stage. 



l + c 



+ 



1 - c 



L(i + x) 2 (i + xT 



(32) 



3.2.2. Unbounded frozen states (Ai > 0) A second type of time-translation invariant 
solutions can be found upon making the ansatz 

Ci,* = 1, G^ = G(t-0 VM', (33) 

but giving up the assumption that X(t) = Ao be constant. Instead, these solutions 
correspond to a fully frozen state, and we find in numerical simulations that 

X(t) = A + Xit (34) 

for a constant Ai > (yet to be determined) and that the qi(t) take the form = t> ji. 
In order to compute the persistent order parameters, it here turns out to be convenient 
to work with the effective process rather than with the closed equations for C and G. 
Upon making the ansatz q(t) = vt for solutions of the single-effective agent process, 
and assuming that 6(t) = 9 is a persistent perturbation one performs an average of the 
effective process (J10J) over the time t (a similar time-averaging procedure for the effective 
process of the conventional batch MG is discussed in [10J ) and obtains 



9 ~ T7TT — V + + V^7- 35 
Ai(l + x) Ai 
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Here r\ = lim^ooT 1 ^2 t<T f](t) is a static Gaussian variable of zero average and with 
variance (jf)^ = (1 + cq)/(1 + x) 2 (with cq = 1 in the frozen state under investigation). 



The only purpose of 9 is to generate response functions, and noting that derivatives 
with respect to 6 can (up to a factor y/a) be replaced by derivatives with respect to rj 
we may set 9 = and write 



Vat] 



1 + — ( — - 
Self-consistency then demands 



(36) 



Aiv^a \dfj / 
Using (|3Tj) these relations are given by 

Xl = 1+ ^- K )T^rr x = Al ( 1+ .( lfe - K ))- (38) 

and the solutions are found to be 

3 /(J! 1 

Ai = — 1 — a H =- + ock, y = . . (39) 

y/2 y/2a-l 



a 2 = - 
2 



In particular this means that the assumption of finite integrated response \ breaks 
down at a = a C 2 = 1/2 so that we expect a transition to a non-ergodic phase. Note that 
this happens irrespective of the value of k, in particular we could not find any signs of 
memory onset at finite x m the case k > 0, as it is observed for the conventional MG 
with impact correction 

The volatility in the frozen ergodic phase is obtained by inserting the above 
expression for \ into (}3*2*|) (with Cq = 1), and one finds 
1 r 1 i 2 

V2~^ ■ (40) 
L y/a\ 

This expression is exact and valid in the interval a G [a C 2 = 1/2, a^n) for any fixed 
value of k. The zeroes of Ai = Ai(a) are found as: 

5 + 4k ± 3 VI + 8k , s 

a = — - 9 . 41 

4(1 - k) 2 

While the upper value agrees with a c \ marking the boundary between the frozen and 
oscillatory phases with finite integrated response, we find that the lower zero of Ai has 
no relevance for k > as it occurs at a lower value of a than the onset of anomalous 
response at a C 2 = 1/2. 

We have tested our theoretical predictions for the order parameters in the ergodic 
states against direct numerical simulations of the batch process. Results for the game 

+ We have performed a calculation along the lines of 0, adapted to the present model, and could not 
observe any continuous onset of long-term memory. Also our simulations show near perfect agreement 
with the ergodic theory for k > so that we have no reason to suspect the breakdown of the assumption 
of weak long-term memory. 




Figure 1. Volatility er 2 , oscillation amplitude c% = 1 — Co, and order parameters 
Ao and Ai for the game without impact correction (k = 0) for different initial 
conditions. Markers are from simulations with go = 0.01 (circles) and qo = 3.0 
(squares) respectively. Solid lines are the theoretical predictions in the ergodic phases 
(continued as a dashed line into the non-ergodic phase for Ao). O labels the oscillatory 
phase of section 13.2.11 F the frozen unbounded phase (section 13.2.2(1 and AFO the 
phase with anomalous response, in which oscillatory and frozen states can be found 
fsection l3.3|) . The vertical dashed lines mark the locations of the 0<->F and F^->AFO 
transitions, a 2 — 1 is the random trading limit. Simulations are for N — 500 agents, 
run for 2000 batch steps (additionally some runs were performed for up to 20000 steps 
to confirm equilibration in the non-ergodic phase). All data are averages over 10 
realisations of the disorder. 

without impact correction (k, = 0) are displayed in Fig. ^ while Fig. [21 is concerned with 
the game with impact correction (k > 0). In both cases we find near perfect agreement 
between the analytical theory and simulations, deviations close to the transitions points 
are due to finite-size and/or finite running time artifacts. The resulting phase diagram 
will be discussed in more detail below. 



Stationary states of a spherical Minority Game with ergodicity breaking 



12 




a \ \ a 

a c i(/c) a c2 

Figure 2. Volatility er 2 , persistent correlation cq, order parameters Ao and Ai for 
the game with impact correction. Results are shown for k = 0.1 (circles), k = 0.25 
(squares), k — 0.5 (diamonds) and k = 1.0 (triangles). Symbols are from simulations 
with N — 500 agents, 10 realisations, initial conditions qo = 0.01. Simulations are run 
for 2000 batch steps (20000 steps for the data for cq at a < 0.5 in order to minimize 
artifacts due to incomplete equilibration). Solid lines: theory. We show the results 
for a 2 and Cq for the game at k — as dot-dashed lines in the upper two panels for 
comparison. 



3.3. The non-ergodic phase 

Throughout the non-ergodic phase one has x — 00 > an d the stationary state depends 
on the initial conditions from which the dynamics is started. Hence, computing the 
asymptotic order parameters would require a full study of the intermediate transient 
dynamics to keep track of the influence of the starting point. This is beyond the scope of 
the present paper; we are here only concerned with the stationary states themselves. In 
numerical simulations of the non-ergodic regime of the game without impact-correction 
these states turn out to be (i) bounded and (ii) either frozen or oscillatory, so that we 
must assume that the ansatz Ai = and C(r) = c + (1 — co)(— l) r is still valid below 
a C 2 = 1/2. Setting % = oo in Eq. (|32j) (which strictly speaking has been derived only 
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Figure 3. a 2 vs Co for the game without impact correction (k = 0) at fixed a = 0.25 
(circles) and a = 0.1 (squares) respectively. Symbols are from simulations, every 
marker corresponds to a different starting point qi(0) — qo for all i with qq varying 
from 0.01 to 2.0 in small steps. Leftmost markers correspond to qo = 0.01, rightmost 
symbols to frozen states resulting from high initial bias. Solid lines are Eq. JJUJ. 



under the assumptions of ergodicity and time-translation invariance) leads to 



For the frozen states (cq = 1) with diverging integrated response we thus expect the 
volatility to vanish as confirmed in numerical simulations, see Fig. ^ Furthermore, one 
realises that in the derivation of (j31j) only the case uo = it in Eqs. ()21J22j) is used, so that 
x' and Ao can be computed independently of c u and x- We may hence conjecture that 
(13 1|) still holds in the non-ergodic regime (at least in approximation). Note however that 
(13 ip can only be expected to be valid in oscillatory states as we make explicit use of the 
assumption c\ > 0. Fig. [T] confirms that A = (a + 1 + 2y / a)/2 appears to describe the 
simulational results also in the non-ergodic regime, as long as the system is initialized 
at a small bias leading to an oscillatory regime. Insertion into ()42|) gives 



As shown in Fig. El this relation turns out to be valid to a very good accuracy. The set 
of possible solutions in the non-ergodic phase thus appears to be parametrized by Co, 
with different values attained depending on initial conditions.* 

* Simulations of the game without impact-correction and at fixed a < a C 2 = 1/2 seem to indicate that 
oscillatory solutions are realised for initial biases qo smaller than some critical value qo c (a). In this 
regime Co is an increasing function of qo- Starting points qo > qoc{oc) lead to frozen states with Co = 1 
and vanishing volatility. From our simulations we cannot rule out that qo c {a) may indeed depend on a. 
A similar threshold separating non-ergodic states with diverging and vanishing volatilities respectively 
has been identified in ^U] in the limit a — > of the conventional MG. 




1 l-co 



(42) 




(43) 
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Figure 4. Phase diagram in (a, «;)-plane. O marks the ergodic oscillatory 
regime at finite integrated response, F is the unbounded, but ergodic frozen phase, 
AF the phase with anomalous response at a < a C 2 = 1/2 and k > in which 
the system is frozen. Finally AFO marks the phase with anomalous response at 
k = 0,a < ac2 = 1/2, in which both frozen and oscillatory solutions are found 
depending on initial conditions. 



For the game with impact-correction (k, > 0) we could not find any oscillatory 
solutions below a = 1/2 in numerical simulations. Independent of initial conditions, all 
observed states were fully frozen and unbounded in the sense of the above definitions, as 
indicated by non-zero values of Ai in Fig. |21 and showed a zero volatility for all a < 1/2. 

4. Phase diagram and discussion 

The resulting phase diagram is depicted in Fig. and we will now turn to a brief 
discussion of its main features, treating the two cases k = and k > separately: 



For the game without impact-correction, we find three distinct phases: (i) an ergodic 
phase with oscillatory behaviour of the correlation function at a > a c \ = 2 (O), (ii) an 
unbounded frozen (but still ergodic) phase at intermediate a £ [1/2,2] (F), and (iii) 
a phase with anomalous response below a = a C 2 = 1/2 , in which we can find both 
oscillatory and frozen states (AFO) , depending on initial conditions (as illustrated in 
Fig. [TJ. Note also that the volatility displays a characteristic minimum at a = a c = 1/2, 
very much like the fluctuations of the original MG (we find a 2 (a c ) = 0, however, whereas 
the volatility in the conventional MG is always strictly positive at finite values of a). 
In contrast to the different spherical model discussed in P we find no discontinuities of 
the volatility or oscillation amplitude across any of the transitions. 



4.1. k = o 
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4.2. K>0 

The behaviour of the model with impact correction is qualitatively similar to the case 
k = for large values of a and as long as k 7^ 1. There we find an oscillatory phase 
at a > a c i(fv) separated from a frozen (but ergodic) phase at intermediate values 
of 1/2 < a < a c i(n). The critical value a c i(n) increases as k is increased, with 
lim K _>i a c i(n) = 00, see Eq. (]2*%|) . so that no oscillatory phase is found for the game 
with full impact-correction (/c = 1). Simulations confirm indeed that the game with full 
impact correction is in fully frozen state (cq = 1) for all a , see Fig. in analogy with 
the conventional MG, where the same feature is observed at k = 1 [8 . We note however, 
that the spherical model does not appear to exhibit a simultaneous onset of memory 
and replica symmetry breaking at finite susceptibility x, as found in the conventional 
MG for < k < 1 pj. Instead we find that anomalous response (x — > 00) sets in at 
a = 1/2 irrespective of k. 

Below a = 1/2 the system is hence in a state of anomalous response both for k = 
and for k > 0. However, we do not observe any oscillatory states in simulations in this 
regime for k > 0, while both oscillatory and frozen states are found below a = 1/2 in 
the game without impact correction. In the above phase diagram we refer to the phase 
with anomalous response as AF for k > 0, as opposed to AFO for k = 0. In particular 
the volatility a 2 turns out to vanish in this AF phase at k — 0, whereas oscillatory 
high- volatility solutions can be found for k = 0,a < 1/2 and small initial bias. This 
results in a discontinuity of a 2 at k = + at any fixed value a < 1/2 (and for sufficiently 
small initial bias), as shown in the upper panel of Fig. El We here plot a 2 as a function 
of k at fixed a = 0.1 and find that o 2 = for all k > 0, but that a 2 is positive at k = 0. 
The magnitude of this discontinuity depends on a and increases as a — > due to the 
increasing volatility in the oscillatory states of the game without impact-correction in 
this limit (Fig. Note that the volatility is continuous at k = + in the ergodic regime 
a > 1/2, see the lower panel of Fig. Again very similar behaviour has been found in 
the conventional MG and also in MGs with dilution [TH ITT]. 

A further parallel to conventional MGs can be drawn upon noting that the 
normalization parameter X(t) = X + Ait remains finite asymptotically (Ai = 0) for all 
< k < 1 above a c \{K) but becomes divergent (Ai > 0) in the F phase at intermediate 
a. This reflects the behaviour of the fraction of frozen agents in the conventional MG 
(corresponding to runaway solutions) which for k = is known to rise sharply as a 
is increased from the non-ergodic to the ergodic region. It decreases as a is increased 
further for conventional MGs at all < k < 1 and vanishes in the limit a — > oojj. In 
addition the fraction of frozen agents approaches unity as a — > for the conventional 
MG with impact correction (« > 0), in analogy with the observed non-zero values for 
Ai in the AF phase of our model with n > 0, see Fig. 121 The spherical model with full 
impact correction (k = 1) has Ai > for all a, its conventional counterpart is in a fully 

(t Note however that simulations of the game at k = in the AFO phase indicate that Ai = even for 
high bias. In the conventional game one here hnds a large fraction of frozen agents. 
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Figure 5. a vs k for different values of a at initial bias qo — 0.01. Upper panel: 
a = 0.1 (circles), where the system is in a non-ergodic state for all k. Lower panel: 
a = 3.16 (diamonds) and a = 1 (squares). For a = 3.16 the system is in the F phase 
above a critical value of k « 0.12, and in an oscillatory state O below n w 0.12. For 
a = 1 the F is realised at all n. Both panels: markers are from simulations, solid lines 
from theory for the respective phases. The dashed line in the upper panel is a guide 
to the eye. 



frozen state irrespective of a. 
5. Replica theory 

We will now turn to the statics of the model. The starting point of the analysis of the 
statics of the standard MG without spherical constraints is the random function 

H «<V) = ^N^ (E^V^ + Sv^ + NWW + (1 - ^ ) '( 44 ) 

\x \ i,j i i J 

The computation in the standard MG then proceeds by a minimization of H K in terms 
of the variables {tpi} using replica techniques. The {y?i} here correspond to temporal 
averages over the (discrete) spins Si(t) e { — 1, 1}, i-e. one has tpi = lim^oo t~ l ^ T<t Sj(r) 
so that the {tpi} are continuous variables —1 < ipi < 1. In the spherical model 
we have Yli 0«(^) 2 — N at any given time t so that the temporal averages y?j = 
lim t ^oo t^ 1 Yl,T<t ( Pi( T ) are no ^ individually constrained to the interval [—1,1], but fulfill 
a global constraint of the form £V </?f < iV. That is to say the vector cp = (</?i, ...,ip N ) 
can take values corresponding to any point within a sphere of radius v^/V, and is not 
constrained to its surface. Note that this inequality holds as an equality only if all 
{4>i{t)} are time-independent, i.e. if a fully frozen state is reached. To capture cases 
different from this fully frozen state we first constrain <p to a shell of radius ry/N, and 
will subsequently minimize the free energy with respect to < r < 1. We expect to 
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recover the frozen states as minima of H at r = 1, and will see below that minima at 
values r < 1 correspond to the oscillatory stationary states of the batch game. 

At fixed r the replicated partition function at an 'annealing temperature' T = 
is given by 



{Z r 



r 2 N 



^T(y^) 2 exp [-(3H K (^)} 



(45) 



where we have introduced replicated fields a — 1, . . . ,n is a replica index, and 

ip a = (</?",..., tp%). The delta-functions in the above partition function enforce the 
spherical constraints on each of the sets {v^}, a = 1, . . . ,n. 

We will be interested in the limit (3 — > oo of the above expression eventually, as the 
minima of H K are found as 



lim min 



lw> — — = ~ nm li m li m — • (46) 

N^oo x ^ s N ^ooJV^oon^O (3Nn 

The notation H K (ip)\r) indicates that the fields are restricted to radius r at this stage. 
The computation of (Z r ) n is lengthy but standard and follows the lines of [T3J[TT]. Upon 
making a replica symmetric ansatz one ends up with the following expression for the 
free energy given r 



Z n - 1 



fr,Rs(Q, q, S, S, p) 



a 
2P 



log 



0, 



i + HQ-q) 

a 



+ 



a 



2a + f3(Q-q) 



sq) 



1 



K , 



Q) + 



log / d(pexp(-/3V z (ip)) 

Q), 



(47) 



where Q and q are the diagonal and off-diagonal elements of the replica symmetric 
overlap matrix respectively, and S and s their conjugate variables, p is a Lagrange 
parameter enforcing the spherical constraint Q = r 2 . V z (ip) is an effective potential 
defined by 

V z (cp) = -V^sz V - ^(S - s)v 2 (48) 

and (. . .) z denotes an average over the standard Gaussian variable z. 

The corresponding saddle-point equations are obtained by working out the 
variations of f r ,RS with respect to the parameters {Q, q, S, s, p}. By construction, 
df r ,Rs/dp = implies Q = r 2 . The other saddle-point equations are found along 
the lines of ^3 [n]. In the limit (3 — » oo we will be looking for solutions with 
linig^oo q = lim^^oo Q and lim^oo s = lim^oo S, and one defines the finite quantity 
X = lim^oo -{Q — q). After some algebra the saddle-point equations lead to 

X =[^l + l/r*-l}-\ (49) 

in direct correspondence to (J27|) . The free energy at zero temperature reduces to 



lim f r RS 

fj— >oo 



1 1 + r 2 

2(I+xF 



K , 

+ 2 (1 



r 



Vl + r 2 — r/y/a 



2 + ^(l-r 2 ). 



(50) 
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Finally we minimize this expression with respect to < r < 1. If the minimum is taken 
at r = 1 we expect the system to be fully frozen, while a minimum at r < 1 corresponds 
to an oscillatory state as found in the dynamics. For a > a c \ (with a c i as given in Eq. 
()28))) the absolute minimum of (|5U|) in the interval r e [0, 1] is identified as 

2 _ 1 + (1 - n)a + a/P + (1 - K ) a ? ~ 4a 

For a < a c i the minimum is attained at r = 1. Hence replica theory leads to the same 
expressions for the persistent order parameters as the generating functional analysis (see 
e.g. Eq. (f2"Kf )■ and in particular we find the same phase boundary between the O and 
F phasesff. Furthermore, the frozen phase (r = 1) obtained from the statics breaks 
down when x — ¥ i- e - at a = a C 2 = 1/2 (see Eq. (jjSj) ), again in agreement with the 
phase diagram obtained from the dynamics. 



< 1. (51) 



6. Continuous-time limit and general time-step 

6.1. Continuous time 

In this section we will consider a continuous time limit of our model. For simplicity we 
will restrict the discussion to the case without market impact correction, a generalization 
to non-zero values of k is straightforward. We will thus consider 

^qM=-hi-Y ( Jij<j>M ( 52 ) 

3 

where &(t) = ft(t)/A(t), with A(t) > chosen such that £*Li 0< (*) = N ++ - The 
effective single-agent process corresponding to reads 



q(t) = ~aj dt'(l + G)-\t, t')<p(t') + 9{t) + Va7j(t) 



(53) 



(see also jTHI for p-spin spherical spin glass models in continuous time). The order 
parameters C, G and A are to be determined self-consistently according to 

(v(tW)l = [(1 + G)- X D{1 + G T y 1 ] (t, 0, (54) 

C7(t, = mm'))* , G(t, if) = (4>{t)) M , C(t, t) = 1. (55) 

Conversion into closed equations for the correlation and response functions gives 

^L[\(t)c(t,t')} = -^((i + Gr^M') 

+ a((l + G)" 1 D(l + G' T )- 1 G' T ) (t,t'), (56) 

ff Note that while the static approach allows one to distinguish between frozen phases (cq — 1) and 
phases which are not fully frozen (co < 1), the shape of the correlation function cannot be determined 
from the replica analysis presented here. The referral to the latter phase as 'oscillatory' is therefore 
motivated by the previous analysis of the dynamics of the batch game 

++ Notc that this is similar but not identical to a Langevin equation qi(t) = — with H = 



Stationary states of a spherical Minority Game with ergodicity breaking 19 

^ [X(t)G(t, t')} = - a ((1 + Gy'G) (t, t>) + S(t - t'). (57) 

Assuming that a time-translation invariant bounded state is reached, we may again write 
these equations in terms of the Fourier transforms of C and G, and find (A(t) = Ao) 

%u\qC(uj) = - a \-a v '— v 1 (58) 

1 + G{u) |1 + GM| 2 

iuX G(u) = - a-^j- + 1 (59) 
1 + G{uj) 

Using D(t — t') = l + C(t — t') and writing \ — G(0) as usual, we find that the equation 
for C gives 

(zw\l + G(uj)\ 2 + aj C{uj) = a X 27T5(uj). (60) 

Similarly to calculation in the spherical MG in discrete time, we conclude that C(u) 
can be non-zero only for u = 0. In contrast to the model in discrete time, there are 
no oscillatory solutions. This does not seem to be surprising since the continuous- 
time dynamics of (|52)i can be expected to effectively correspond to a system with 
asynchronous random updating (as opposed to the parallel dynamics considered earlier) 
for which also in the conventional MG model global oscillations are suppressed j!6j . 
While this correspondence to asynchronous updating in the MG still deserves further 
attention we note that synchronous updating in neural networks results in a description 
of the effective dynamics in terms of a discrete-time process, while asynchronous 
updating leads to equations in continuous time; see for example [13 CHj for further 
details. 

It is easy to check that assuming a fully frozen state (co = 1) in Eqs. (|HHj) and 
(fSTIj) (evaluated at u = 0) implies a = 2, so that we conclude that there can be no 
extended phase with bounded values of X(t) in the model in continuous time, neither a 
frozen nor an oscillatory one. It would then appear sensible to revert to our unbounded 
state \(t) = A + Ait, with Ai > and to assume that <p(t) in the effective process 
approaches a fixed point (with value dependent on the realization of the single-particle 
noise). This then leads to Eq. (|35|) (where now v = lim^oo q), and hence to exactly the 
same solutions for % and a 2 as in the F phase of the game in discrete time. In particular 
one expects a breakdown of ergodicity at a = 1/2 also for the model in continuous time. 

6.2. Arbitrary, but finite time-step 

One may ask whether this model in continuous time can be obtained as the limit 5^0 
of a model with time-step 5: 

* (t+ *|-* (t) = -ft,-E-w>- («) 

j 

5 can here also be understood as a learning rate, with which the agents process 
information and adjust their score updates and trading actions. We note that the left- 
hand side can be written as [\(t + 8)4>i(t + 5) —\(t)</>i(t)]/5 upon using (f>i(t) = qi(t)/X(t), 
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so that the only effect of the arbitrary time- step is a re-scaling of the {A(t)}. Numerical 
simulations (not shown here) confirm that the oscillation amplitude of the correlation 
function C(t — t') = (<p(t)(j)(t'))^ in the stationary state and the volatility are not affected 
by the introduction of a general time-step, but that the only effect is a rescaling of the 
{<&(£)} relative to the {<f>i(t)}. This is at variance with the spherical model of [T], in 
which the phase diagram changes as the time-step (or equivalently the learning rate) is 
modified, as sketched in the appendix. Note that a subtle dependence of fluctuations 
on the learning rate has been observed in a different agent-based model with continuous 
degrees of freedom in [T5J |20] . 

7. On-line game 

Finally, it is interesting to compare the stationary states of the on-line model (defined 
by Eq. (JJJ)) with those of the batch game. The persistent order parameters (such as 
Co and x) are known to be identical in the stationary ergodic states of conventional 
batch and on-line MGs respectively [Hj, but differences in the volatility are in principle 
to be expected, as the latter depends not only on the persistent parts of the dynamical 
order parameters, but also explicitly on their transients. While these differences between 
batch and on-line games are negligible in conventional MGs with uncorrelated strategy 
assignments, crucial differences in their quantitative values and qualitative behaviour 
were reported for anti-correlated strategies in ^H]- These differences can be attributed 
to persistent oscillatory behaviour of the batch game. 

We compare results from numerical simulations of the batch and on-line games in 
Fig. |H1 We restrict the discussion to the case without market impact correction here 
(k — 0). As expected we observe an ergodic/non-ergodic transition at a = a C 2 = 1/2 
also in the on-line game and find that the persistent correlation Co of the on-line game 
takes the same values as that of the batch game in the ergodic regime. In particular 
Co = 1 in both games in the fully frozen state at intermediate a. Within the accuracy of 
the simulations the locations of the transitions are insensitive to the timing of adaptation 
(on-line versus batch) and in agreement with the batch theory. 

The volatility cr 2 , however, shows a more intricate behaviour, see right panel of 
Fig. El While the volatilities of the two games agree perfectly in the frozen ergodic 
phase at intermediate a G [1/2, 2] and are there given by Eq. PUJ) . one finds differences 
between the batch and on-lines cases in the regime of high values of a, where the batch 
system oscillates in its stationary state. In this phase at a > 2 the volatility of the 
batch game is given by Eq. (|H2*|) . and contains contributions from the persistent as well 
as the oscillatory parts of the correlation and response functions. As shown in Fig. El 
the volatility of the on-line game is consistently smaller than that of the batch game 
in this regime. The spherical game presented here therefore presents another example 
of a MG in which the market fluctuations of the batch and on-line games differ from 
each other. As in these differences can be attributed to the presence of persistent 
oscillations present in the batch game, but not in its on-line counterpart. Indeed, the 
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Figure 6. Persistent correlation Co and volatility a 2 for the batch and on-line games 
without impact correction [k — 0). Open circles are obtained from simulations of the 
batch game (simulation parameters as in the other figures) , open squares correspond to 
the on-line game (simulations performed at aN 2 = 10 4 , run for 100000 on-line steps, 
20 samples of the disorder), qo = 0.01 in both cases. Solid squares correspond to the 
on-line game with initial bias qo = 3. Solid lines are theoretical predictions, with the 
upper line for a 1 in the high-a phase corresponding to the batch theory (Eq. (|32() ') 
and the lower solid line to the on-line approximation l|62|) . The expression for the fully 
frozen phase has been continued into the high-a regime (where it is no longer valid 
neither for the on-line nor the batch game) as a dashed line. 



volatility of the on-line game is consistently larger than the theoretical prediction one 
would obtain in a fully frozen state, see Fig. El and the now standard approximation 

1 + co 



2 



L(l + X) 



+ (1-co 



(62) 



for the volatility in on-line games, obtained by suitably removing all non-persistent 
contributions to C and G (for details see [ZD El) gives near perfect agreement with the 
numerical data. While ()62|) is at present only an approximate expression, an exact 
dynamical theory appears feasible for the present spherical game also with on-line 
dynamics (for example along the lines of [2*T] ) and might allow one to derive an exact 
result for the volatility. Finally, we note that started from equal initial conditions the 
volatility of the on-line game is consistently smaller than that of the batch game also 
for a < 1/2 (if the initial bias in this non-ergodic regime is small enough not to end 
up in a fully frozen regime, where a 2 = for both games). We recall that under these 
circumstances the stationary state of the batch game is again an oscillatory one. 



8. Concluding remarks 

In summary we have presented a detailed analysis of the stationary states of a spherical 
version of the MG, both from a static and a dynamical point of view. While the present 
model defines only a minor modification of an earlier spherical MG, we find that some of 
its features bear more similarity with those of the original MG. In particular one finds 
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a non-ergodic regime, characterized by a strong dependence of the stationary states 
on initial conditions, while the spherical model of pQ turned out to be ergodic for all 
values of the model parameters. Also the model presented here does not display any 
discontinuities in the volatility, but the characteristic minimum of a 2 at the point at 
which the onset of non-ergodicity occurs. Despite this increased resemblance to the 
original MG, the model is still exactly solvable in the ergodic phases, and in particular 
exact expressions for the volatility can be obtained. The model also allows to study 
the dependence of volatilities on the timing of adaptation of the agents. Similar to [TH] 
one finds that batch games in their oscillatory phases exhibit larger fluctuations of the 
total bid than an on-line game with the same model parameters. We are here able to 
give an exact expression for the batch volatility, and a very accurate approximation for 
the on-line case, whereas in [TH] analytical results for the volatility of the batch game 
with fully anti-correlated strategies were restricted to upper and lower bounds. Hence, 
the present spherical model might serve as a starting point for a more detailed study 
of the differences between on-line and batch games, which are reminiscent of the ones 
recently observed for example in [T%1 12*2~] for sequential and synchronous updating in 
other spin systems. Two-cycles, i.e. persistent oscillations, play a crucial role in those 
studies similar to what is observed in our model and in the MG with anti-correlated 
strategies. A detailed analysis of an on-line version of the two spherical models might 
also help to elucidate the role of transients in the dynamics and the success of the 
ad-hoc removal of non-persistent parts of the correlation and response functions when 
deriving an approximate expression for the volatility. Finally, the dependence on initial 
conditions observed in the model discussed here might motivate further studies of its 
non-ergodic regime. Ideally, one would like to be able to compute the persistent order 
parameters as a function of the initial bias qo, which would most likely involve a detailed 
computation of the transient behaviour. One might hope that work along these lines 
could also help to shed some light on the non-ergodic, turbulent regime of the original 
MG, a thorough analytical understanding of which is still awaited. 
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Appendix: Comparison of the two spherical models 

In pQ a spherical limit of the dynamics of the MG was taken upon using the update rule 
[1 + X(t + l)] qi (t + 1) = qi (t) -hi Jijqjit), (63) 

3 

where the {A(£)} are chosen to impose a global spherical constraint ^2 i=1 qi(t) 2 = Nr 2 . 
A perturbation field 6{t) may be added to the right-hand-side to generate response 
functions. The control parameters are given by a and r, and in order to be consistent 
with the spherical constraint initial conditions {qi(t = 0)} were chosen such that 
2\^i<?«(0) 2 = Nr 2 . In this appendix we will first briefly summarise the main findings 
of pQ and then extend the analysis of 1J by studying the approach to its continuous- 
time limit and by a brief discussion of the statics of this earlier spherical model. 

Phase diagram 

Similar to the analysis in the main part of the paper the generating functional 
calculation in pQ leads to a coupled set of closed equations for the correlation and 
response functions (note that here one has C t t> = liniAr^oo iV _1 £\ (qi(t)qi(t')) and 
G w = lirn^oo N' 1 djq^tj)/d9(t>)): 

[1 + X(t + 1)] C w = Cte +«[(!+ G)- 1 J D(1 + G T )~ 1 G T ) U > 

-aKl + G)- 1 ^, (64) 
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Figure 7. Phase diagram of the spherical MG studied in 1 , displaying three phases: 
O, oscillating correlation function and finite integrated response; F, frozen phase with 
finite integrated response; and AF, anomalous frozen phase with diverging integrated 
response. Throughout the phase AF the volatility is zero. The 0<->F transition and the 
F^->AF transition are continuous. The discontinuous transition from 0<->AF occurs at 
a = 3-2^2 ~ 0.172. The triple point corresponds to a = 3-2^/2 and r = r* w 0.455. 
The dotted-dashed line separating the two frozen phases is given by a — r 2 / (1 + r 2 ) . 



[1 + X(t + 1)] G w , = G u > - a[(I + G) G]tf + 6 W , (65) 

from which the stationary states of the model can be computed in a very similar way 
to the analysis in the main text of this paper. The resulting phase diagram is depicted 
in Fig. [3 and displays three distinct phases O, F and AF, similar to the model studied 
in this paper. All phases are bounded however, as runaway solutions \qi(t)\ — > oo 
are impeded by the spherical constraint on the {qi}. While the transitions 0^->F and 
F^AF are continuous, a jump occurs in the volatility and oscillation amplitude across 
the O^AF transition. In particular one finds that the volatility vanishes identically 
throughout the AF phase, where x — oo. The macroscopic order parameters (including 
a 2 ) show no dependence on initial conditions in any part of the phase diagram, in 
contrast to the model devised in this paper**. 

**Note that in [Q the second moment of the distribution of initial conditions was fixed by 
iV _1 X^*(0) 2 = r 2 for consistency as this normalization on the {<&(£)} is imposed for all subsequent 
times. The freedom of choosing initial conditions in £Q was thus limited to distributions with this fixed 
value of their second moment which in turn plays the role of a model parameter. In the present model 
no constraint on the {<&(£)} applies so that initial conditions {qi(t = 0)} can be chosen arbitrarily. The 
analysis in the main part of the paper shows that different second moments A(0) 2 = TV -1 J2i 9i(0) 2 
lead to different macroscopic stationary states for a < 1/2. 



Stationary states of a spherical Minority Game with ergodicity breaking 
4 



25 




















7 F 

















a 



Figure 8. Phase diagram of the model of for general time-step 5. Left (5 > 1/2): 
solid lines: continuous O^F transition for 6 = 1.0,0.8,0.6,0.55 (from left to right). 
The dotted-dashed line is given by a — r 2 /(l + r 2 ) and marks the F^AF transition. 
Right: phase diagram for 5 < 1/2: solid lines: continuous O^F transition for 
5 = 0.4, 0.3, 0.2 (from left to right). The dotted-dashed line is given by a = r 2 /(l + r 2 ) 
and marks the F^AF transition. 



Approach of the continuous-time limit 

The model defined by (|63p is easily generalised to arbitrary (but finite) values S of the 
batch time-step: 

[i + A(«+')y+*)-*m = _fr_g jBaW , m 
^E«. 2 w = r2 fora1 "- < 67 ' 

i 

and the analysis of the stationary states is easily adapted accordingly. In contrast to the 
model discussed in the main part of the paper, one here observes interesting effects on the 
phase diagram. For l/2<5<2we find that the phase diagram and types of transitions 

are essentially preserved, but that the triple point moves according to a = (^\/2/5 — lj 
along the line a(r) = r 2 /(l + r 2 ) as 5 is varied, and that the numerical values of the 
locations of the different transitions change accordingly, as depicted in the left panel of 
Fig. |HJ In the case of small 5 < 1/2 the discontinuous 0<->AF transition is no longer 
present, and the range of the oscillatory phase moves towards large values of a and r 
as S is decreased, see the right panel of Fig. |HJ In the limit <5 — ► the oscillatory phase 
will eventually no longer be present, so that the resulting model in continuous time only 
displays the F and AF phases, respectively, separated by a continuous transition at 
a = r 2 / (1 + r 2 ). As before all order parameters can be computed exactly for arbitrary 
values of 5 in all phases, and numerical simulations (not shown here) confirm these 
results as well as the changes in the phase diagram 



Statics Finally, we note that the statics of the model (jo^j) can be studied upon 
minimizing the Hamiltonian defined in 1)44)1 2TH. Again the variables {<^} are 

to be constrained to a sphere of radius r, where now r is a model parameter and 
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fixed (i.e. the free energy is not minimized with respect to r). Indeed the dynamical 
approach via generating functionals and the static calculation via replica methods turn 
out to deliver the same order parameters in the F phase of the model defined by (JfiHjl . 
However, a replica approach along the lines of the main part of this paper (in which 
ground state solutions are not constrained to the surface of the sphere with radius r 
but can lie anywhere within) fails to reproduce the dynamical results obtained for the 
oscillatory regime of ([63)1 . It appears that here a more sophisticated approach based on 
(pseudo-) Hamiltonians of the Peretto type (2H ESj might be required, possibly along 
the lines of [HI 122- 



